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SUMMARY 

Currently, the heat transfer equation used in the rotary combustion engine 
(RCE) simulation model is taken from piston engine studies. These relations 
have been empirically developed by the experimental input coming from piston 
engines whose geometry differs considerably from that of the RCE. The objec- 
tive of this work was to derive equations to estimate heat transfer coeffici- 
ents In the combustion chamber of a RCE. This was accomplished by making 
detailed temperature and pressure measurements in a direct injection stratified 
charge (DISC) RCE under a range of conditions. For each specific measurement 
point, the local gas velocity was assumed equal to the local rotor tip speed. 
Local physical properties of the fluids were then calculated. Two types of 
correlation equations were derived and are described in this paper. The first 
correlation expresses the Nusselt number as a function of the Prandtl number, 
Reynolds number, and characteristic temperature ratio; the second correlation 
expresses the forced convection heat transfer coefficient as a function of 
fluid temperature, pressure, and velocity. 


INTRODUCTION 

The results of heat transfer studies over the past several decades 
(refs. 1 to 5) have been applied to internal combustion engine analysis to 
access engine performance and structural integrity. Within the past 5 yr, 
extensive implementation of this work in computer codes has been used to simu- 
late piston engine operation. In general, the equations used for heat transfer 
analysis in the piston engine have taken the following form: 
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where equation (1) is used to estimate the heat transfer coefficient in the 
reactor (ref. 4); equation (2) is used to estimate the turbulent heat transfer 
coefficient through a flat plate (ref. 8); equation (3) is the equation derived 
by Woschni (ref. 5). 

These equations have been found to be satisfactory when the constants in 
them were based on empirical information available for the piston engine of 
interest. 

Previous simulations (refs. 7 to 9) of the heat transfer processes in the 
RCE have assumed that the turbulent gas-to-wall heat transfer coefficient is 
given by the McAdams correlation for turbulent flows in pipes; however, these 
assumptions were made without the benefit of the availability of detailed 
experimental data. Data to make this correlation possible was collected by 
NASA during this analysis. 

The applicability of the above heat transfer equations to the rotary com- 
bustion engine is described in this study. 


ROTARY ENGINE TEST RIG 

The test rig (fig. 1) was constructed using a single rotor, rotary engine 
built by OMC, with a 39.6 in. 3 displacement. The engine is naturally aspirated 
and incorporates direct fuel injection combined with spark ignition to create a 
stratified charge environment for combustion. 

One fuel injector and one spark plug were used in this configuration. The 
rotor was cooled by the Intake air seeded with oil to eliminate the need for a 
main bearing lubrication system. End bearings were lubricated in a conven- 
tional fashion. Coolant is directed peripherally around the hot side of the 
rotor housing to simplify the rotor housing and the side housing geometry. 
Figure 2 shows the path of the coolant through the engine. 

The fuel system consisted of a distributor type injection pump and a four- 
hole pencil injector which was located near the minor axis of the trochoid 
housing. The fuel used was an ASTM diesel check fuel with a cetane number of 
53.7. 


A single spark plug located on the minor axis of the trochoid housing was 
used to initiate combustion. The single spark per power stroke was generated 
by a capacitive discharge ignition system. 

Lubrication for the main bearings, rotor bearings, and the rotor seals was 
provided by two solenoid-operated metering pumps. The main bearing lubrication 
oil was a synthetic, two-cycle oil which was not recirculated but consumed in 
the combustion process. The heat of combustion of the oil was negligible in 
comparison to that of the fuel. 

The remainder of the test rig consisted primarily of an eddy current dyna- 
mometer and an electric motor. The motor was capable of motoring the engine up 
to speeds of 3200 rpm, while the dynamometer could absorb engine torques up to 
50 ft-lbf at speeds up to 6000 rpm. 
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INSTRUMENTATION 


The flow of air into the engine was measured using a laminar flow element. 
The flow element was separated from the engine by a plenum chamber which damp- 
ened pulsations from the intake strokes of the engine. 

Fuel flow was measured by a Flo-tron positive displacement flow metering 
system. 

The engine coolant was a one-to-one solution by volume, of ethylene glycol 
and water. The coolant flow was measured by a turbine flowmeter. Both the 
inlet and outlet coolant temperatures were monitored along with the engine 
speed, exhaust temperature and housing temperatures to establish steady-state 
conditions . 

The local average engine temperatures (fig. 3) were determined by the use 
of copper-constantan thermocouples mounted in the rotor housing. A total of 
16 thermocouples were mounted in pairs in eight locations around the trochoid 
housing. At each location, one thermocouple was mounted directly behind the 
other in a straight line which was perpendicular to a tangent to the inner tro- 
choid surface at that location (fig. 4). The second thermocouple was placed 
0.125 in. from the first thermocouple. The thermocouples were positioned in 
this manner so that average local heat flux (fig. 5) and average gas-side wall 
temperature calculations could be made. In addition, all of these thermo- 
couple junctions were located in the midplane of the trochoid housing. Four 
additional thermocouples were mounted in pairs, in two locations in the coolant 
passage, to measure the coolant side surface temperature and the bulk coolant 
temperature. 

Pressures in the combustion chamber of the engine were measured by a 
series of four water-cooled, piezoelectric pressure transducers. These trans- 
ducers were mounted in the midplane of the trochoid housing in such a way that 
the chamber pressure would constantly be monitored (fig. 6). This constant 
pressure monitoring throughout all four engine cycles was made possible through 
the use of a NASA-developed electronic switching and correlating system called 
the Modular Engine Instrument System (MEIS). The composite pressure curve, 
which is a real-time measurement, was recorded on a digital waveform recorder. 
The waveform recorder incorporated a series of individual analog-to-digi tal 
(A-D) converters which permitted a number of individual signals to be recorded 
simultaneously. The commercially available waveform recorder (fig. 7) was mod- 
ified such that it could be triggered by an external signal and then driven by 
an external clock. The recorder was triggered manually when a steady-state 
condition was reached with the engine. The recording was done in a crank angle 
reference frame as opposed to a time-base reference clock. The external clock 
signal and trigger were produced by a NASA designed MEIS/waveform recorder 
interface module. 

Each A-D converter of the waveform recorder had 32K bytes of memory asso- 
ciated with it, which translates into 32 complete, four-cycle, composite pres- 
sure curves. For each data point taken, 32 consecutive composite curves were 
recorded and automatically downloaded to a microcomputer having hard-disk 
memory capability. These 32 composite curves were ensemble-averaged in post- 
run time to form one average composite pressure curve representative of the 
engine conditions at the time data was taken. 
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Exhaust emissions were measured by a Scott Research Laboratories Test Con- 
sole. Sampling was taken at a remote module and pumped 30 ft through a heated 
line to the emissions analyzers. The following equipment was used: 

Scott Model 415 - Total HC Analyzer 
Scott Model 250 - Paramagnetic O 2 
Detector 

Scott Model 125 - NO/NOx Analyzer 
Beckman Model 864 - CO 2 Analyzer 
Beckman Model 865 - CO Analyzer 


METHOD OF CALCULATION 

To simplify calculation, the following assumptions were made: 

(1) Only convective heat transfer was considered. 

(2) Assume no exhaust gas recirculation or inter-chamber gas leakage. 

(3) Uniform chemical composition and gas temperature throughout the cham- 
ber at each instant of time. 

(4) The gas and coolant physical properties can be computed from the aver- 
age film temperature, local pressure, and composition. 

(5) The local gas velocity is assumed to be the average surface tangential 
velocity (it is a time average and is not a spatial average) of the rotor. 

The OMC test engine was run at four selected speeds: 2200, 2500, 2800, 

and 3200 rpm. Temperature and pressure data were taken at 16 thermocouples 
and four pressure transducers (fig. 4) under steady-state operating conditions. 
The local gas composition was obtained by balancing equation (4). Knowing the 
temperature, pressure and gas composition, the local gas properties were 
calculated. 

Gas side convective heat transfer calculations were based on a number of 
empirical and derived equations. Working fluid pressure, exhaust composition, 
and temperature measurements were collected from three different data acquisi- 
tion systems. First, the pressure transducers and MEIS measure and correlate 
local gas pressures. For each crank angle, the pressure is assumed to be con- 
stant throughout the combustion chamber. Second, the Gould 9000 System pro- 
vides the high speed digitization necessary to record the pressure data. In 
the intake and compression processes of the direct injection RCE, the working 
fluid is simply air and water vapor. During combustion, information is based 
on sampling emissions and calculations using the equation: 

^1 ^9 ^19 ^ n 2 ^2 ^ n 3 ^2 n^ H2O “* n^ CO2 + n^ H2O + n^ N2 + n^ O2 

+ n g CO + n ]0 NO + C g H ]g (4) 

The summation of n 5 through nn is the total moles/cycle (Enj) of the 
gas mixture during combustion and exhaust. Third, the values ng through nn 
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are measured by the Scott Emission Analyzer and recorded in the Escort system. 
The local volume of the working chamber is calculated from equation (10): 

_ — 3 VT 2 

V(a) = V m1n + eRb (1 - cos | a) (5) 

The local working fluid temperature can be calculated by the modified 
ideal gas law equation: 


P V(a) = Z ( En i ) R Tg(a) 


( 6 ) 


Using the area surrounding thermocouples 3 and 4 as an example (fig. 4), 
the local average heat loss can be calculated by the equation: 



(7) 


Local gas side average wall temperature can be extrapolated 
from the equation: 



( 8 ) 


Finally, the local gas side heat transfer coefficient is calcu- 
lated from the equation: 



= h (a) (T (a) - T (a) ) 
g g wg 


(9) 


METHOD OF CORRELATION 

The transfer heat in gases has been widely studied (refs. 2 to 5). The 
results of this work have been correlated by a variety of methods which, in 
general, can be represented by two different kinds of correlation. The Corre- 
lation I is to correlate the Nusselt number as a function of the Reynolds 
number, Prandtl Number, and characteristic temperature ratio. Correlation II 
is to correlate the gas side transfer coefficient as a function of gas phase 
temperature, pressure, and velocity. 

The gas velocity v g in equation (10) is made up of three components. 

The first component is caused by the motion of the moving rotor relative to 
the stationary housing (i.e., the combustion chamber is rotating). The second 
component is squishing velocity caused by the chamber's geometry change. The 
third component is the velocity induced by the expanding combustion gases. 

The magnitude and direction of this gas velocity is not known; currently, only 
laser diagnostic techniques can provide this information. Since this informa- 
tion was not available for this study, the gas velocity is assumed to be equal 
to the average (time average at the same crank angle, and is not a spatial 
average) surface tangential velocity of the rotor. This velocity can then be 
calculated by the following equations: 
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V (a) = | co D(a) 

D(a) = 9e 2 + (R 1 ) 2 + 6 eR‘ cos | aj ( 11 ) 

The physical properties used in equation (11) are found for each of the 
eight measurement points along the rotor housing. The locally measured gas 
pressure and housing temperature are used to generate gas temperature (Tg) and 
average gas side wall temperature (T wg ) by equations ( 6 ) and ( 8 ), respectively. 
The local gas composition can be balanced by equation (4), and the local film 
temperature calculated from equation (15). The characteristic diameter (D) can 
be found from equation ( 11 ) yielding enough information to calculate velocity 
(Vg) from equation (10). Thus, at each point all the properties can be esti- 
mated and the gas heat transfer coefficient (hg) from the left hand side of 
equation (14) can be calculated from equation (9). 

The constant terms of the above two correlated equations were derived from 
a multiple linear regression technique. 

The multiple regression determines mathematically how variable y depends 
on independent variables x^ , X£, X 3 , etc. Using matrix mathematics, the rela- 
tionship between the variables is assumed to be linear, and of the form: 

y = a Q + a 1 x 1 + a 2 x 2 + a 3 x 3 ( 12 ) 

The numerical value of the A n are constants. This method provides a 
least-squares best fit by determining the a's for which the sum of the errors 
(i.e., the difference between the actual and predicted y's) is zero and the 
sum of the squared errors is minimum. A coefficient of determination is used 
to indicate how well the entire equation fits, a value of ± 1.0 infers a perfect 
fit. A Lotus 1-2-3 computer program was used to perform this routine. 

Based on the experimental data, a stepwise multiple linear regression pro- 
gram was used. This program sought the most significant parameter from a list 
of candidates provided, and extracted its effect before seeking the next most 
important parameter. The coefficients were not forced, nor was the order of 
parameter selection. 


Correlation I 

The Nusselt number is described as a function of the Reynolds number, 
Prandtl Number, and characteristic temperature ratio. This relationship can 
be represented by an equation of the form. 
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The physical properties involved are evaluated at the gas side film tem- 
perature, which is defined as: 
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f.g 



(15) 


When our experiment data is closely examined, there is a large temperature and 
pressure difference in RCE between unburned gas and burned gas. In order to 
derive a more accurate heat transfer equation, we decided that for each corre- 
lation, the burned gas case and unburned gas case should be analyzed sepa- 
rately, and discussed below: 


Correlation I-A - The Unburned Gas Case . - The experimental data used in 
this unburned gas mixture correlation were taken from the engine motoring cases 
and intake/compression cycle of the engine firing cases. Two regressed equa- 
tions are: 


-4.657 


N = 0.062 ^ 


(Rr)- 4 - 389 (Re) 0 - 745 


-4.552 


N. = 0.3311 ^ 


(Re) 


0.7387 


(16) 


(17) 


The average error of equations (16) and (17) is 15.25 and 15.26 percent, 
respectively (table I). The coefficient of determination of equations (16) 
and (17) is 0.973 and 0.975, respectively (tables II and III). Table I 
reflects that the Prandtl number changes slightly for the unburned gas case. 
This information, along with the higher confidence in equation (17) indicates 
that the Prandtl number is not an important parameter in the unburned gas heat 
transfer equation. 

Correlation I-B - Burned Gas Case . - The experimental data used in this 
burned gas correlation were taken from the combustion/exhaust cycle of engine 
firing cases. Again, the Prandtl number is not a dominant parameter in the 
burned gas analysis. Thus the regressed equation is: 
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13.865 



0.3906 


(18) 


The average error of equation (18) is 22.73 percent (table IV) and the 
coefficient of determination is 0.788 (table V). Note in equations (17) 
and (18) the Prandtl number is not included because its value change is very 
small based on the experimental data; but, if the Prandtl number does vary 
widely, then it should be included in the correlate equation. 


Correlation II 

The gas heat transfer coefficient as a function of gas temperature, pres- 
sure, and velocity. This relationship can be represented by an equation of 
the form: 
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A ' R 1 C 

h = K'(T) rt (P) D (v) L 


( 19 ) 


Again, as mentioned before, the burned gas case and unburned gas case 
were analyzed separately by the multiple linear repression program. 

Correlation II-A - Unburned Gas Case . - Two equations were regressed and 
compared to study whether the bulk gas temperature (T^) or gas film tempera- 
ture (Tf) is a dominant parameter. The equations are: 

h(a) = 19.402(T b )“°- 538 (P) 0,3593 (v) 1 - 2152 <2 

h(a) = 1 .0808(T f )'°’ 092 (P) 0J064 (v) 1 ,2152 (2 


The average errors of equations (20) and (21), are 18.69 and 19.50 percent 
respectively, and the coefficients of determination are 0.9201 and 0.9008, 
respectively. Thus, we concluded equation (20) is better than equation (21) 
and the result is shown in table VI. 

Correlation II-B - Burned Gas Case . - Two equations are regressed to 
determine if either bulk gas temperature or gas film temperature is an impor- 
tant parameter. 


h(a) = 6.2315(T b )" 0 - 45 (P) 0,0689 (v) 1 - 4327 


h(a) = 2.34.1 4(T f ) _0 * 941 (P)°* 1378 (v) 1 * 2767 


( 22 ) 

(23) 


The average errors of equations (22) and (23) are 17.83 and 20.69 percent, 
respectively; the coefficient of determination is 0.8726 and 0.7059, respec- 
tively. Thus we concluded that equation (22) is a better equation to be used 
in the burned gas region. The result is shown in table VII. 


CONCLUSIONS 


1. The investigation of heat transfer mechanisms in the DISC rotary engine 
led to different regressed equations for convective heat transfer. A separate 
equation should be used depending on whether gas is unburned or burned. 


There are two equations derived for the unburned gas region 
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There are two equations derived for the burned gas region: 
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h(a) = 6.231 5(T b )“°' 45 (P)°* 0689 (v> 1 * 4327 


(18) 

( 22 ) 


2. For better pressure resolution, a 10-bit or higher A-D convertor is 
required in order to reduce the experimental error. In this study, the pres- 
sure data digitized by the eight-bit A-D converter does not have adequate reso- 
lution in the lower pressure range. For example, if the maximum pressure in 
the chamber is 3.5 MPa, then the incremental pressure change is 3.5 MPa/(2 8 /2) 

= 0.027 MPa. At peak pressure the error is 3.91/500 = 0.8 percent, but near 
the intake and exhaust the error is 3.91/14.7 = 26.6 percent. 
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APPENDIX 


A Reynolds number exponent defined in equation (13) 

A' temperature exponent defined in equation (19) 

B Prandtl number exponent defined in equation (13) 

B' pressure exponent defined in equation (19) 

b width of the rotor housing, ft or (m) 

C temperature ratio exponent defined in equation (13) 

C' velocity exponent defined in equation (19) 

Cm mean piston speed in ft/sec. or (m) 

Cp specific heat of gas at constant pressure, Btu/(lb)(°F)(or cal /(g) 
(°K)) 

D characteristic length in the gas chamber, ft or (m) 

d characteristic width of the piston, ft or (m) 

e eccentricity of rotor journal, ft. or (m) 

h gas heat transfer coefficient, Btu/(sec)(ft 2 )(°F) or (Kcal/(sec) (m 2 ) 

( °K) ) 

hg(oc) average local gas heat-transfer coefficient, Btu/(sec)(ft 2 )(°F) or 
(Kcal/(sec) (m 2 ) (°K)) as a function of crank angle 

K constant in heat transfer equation, defined in equation (13) 

K' constant in heat transfer equation, defined equation (19) 

k thermal conductivity of gas, Btu/(sec)(ft)(°F) or (cal/(sec) (m) (°K) 

k] average thermal conductivity between points 3 and 4 (Fig. 3), 

Btu/(sec) (ft)(°F) or (Kcal/(sec) (m) (°K)) 

k£ average thermal conductivity between point 3 and gas side surface 
Btu/(sec)(ft)(°F)or (Kcal/(sec) (m) (°K) 

Nu Nusselt number, hD/k 

ni number of moles of component i, Ib-mole/cycle or (g - mole/cycle) 

Enj total number of moles of burned gas, Ib-mol e/cycle or (g - mole/cycle) 

P pressure of the working fluid, psia 

Pr Prandtl number, Cp p/k 


10 



i 
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i 



R 

R' 


local heat transfer rate per unit area through point 3 and 4, 
Btu/ft 2 -hr or (K Cal /hr - in 2 ) 

gas constant, 10.731 (psi )(ft 3 )/(lb-mole)(°K) or (1.987 cal/g - mole) 
( °K) ) 

the generating radius, ft. or (m) 


Re pvD/p 

S 1 distance between points 3 and 4 (fig. 3), ft. or (m) 

s 2 distance between point 3 and gas side surface (fig. 3), ft. or (m) 

T 3 average housing temperature at point 3, °F or (°K) 

T 4 average housing temperature at point 4, °F or (°K) 

Tb average bulk temperature of fluid, °F or (°K) 

Tfg average local gas side film temperature °F or (°K) 

Tg average bulk temperature of gas °F or (°K) 

T w average wall temperature, °F or (°K) 

T W g average gas side wall temperature, °F or (°K) 

Vg characteristic gas velocity, ft/sec or (m/sec) 

V( a) local volume of the working chamber, ft 3 , or (m 3 ), as a function of 
crank angle # 

V m in minimum volume of the working chamber, ft 3 , or (M 3 ) 
p gas density, lb/ft 3 or (g/cm 3 ) 

to angular speed of output shaft, 1 / ( sec ) 

a angle of output shaft 


\ 
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TABLE I 


0.0620 { Tb/Tf ) ~-4.657 (Pr)~-4.389 (ReTO.7450 

0.3311 { Tb/Tf ) ~-4. 552 (Re)~0.7367 
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Number 

Nuf 

Tb/Tf 

Pr 

Re 

Nucal . 

dNu(%) 

abs() 

Nucal . 

dNu(%) 

abs( ) 

1 

254.83 

1.2567 

0.6871 

24 184 

205.01 

-19.55 

19.55 

198.55 

-22.09 

22.09 

2 

136.06 

1.1263 

.6903 

4 664 

98.17 

-27.85 

27.85 

97.24 

-28.53 

28.53 

3 

274.14 

1.2833 

.6882 

27 874 

205.26 

-25.13 

25.13 

200.40 

-26.90 

26.90 

4 

263.76 

1.2788 

.6909 

31 529 

224.81 

-14.77 

14.77 

222.98 

-15.46 

15.46 

5 

100.23 

1.1626 

.6891 

6 221 

105.76 

5.52 

5.52 

104.06 

3.82 

3.82 

6 

350.6 

1.2734 

.6882 

36 931 

262.43 

-25.15 

25.15 

255.40 

-27.15 

27.15 

7 

69.95 

1.175 

.6887 

3 974 

72.27 

3.32 

3.32 

71.27 

1.89 

1 .89 

8 

198.96 

1.3057 

.6873 

23 714 

168.85 

-15.13 

15.13 

164.42 

-17.36 

17.36 

9 

166.53 

1.1126 

.6901 

6 059 

126.46 

-24.06 

24.06 

124.67 

-25.14 

25.14 

10 

297.86 

1.2945 

.6882 

36 017 

238.60 

-19.90 

19.90 

232.66 

-21.89 

21.89 

11 

87.99 

1.1684 

.6889 

5 269 

91 .43 

3.91 

3.91 

90.01 

2.30 

2.30 

12 

255.28 

1.2969 

.6879 

31 391 

213.93 

-16.20 

16.20 

208.48 

-18.33 

18.33 

13 

946.9 

1.172 

.6949 

121 745 

900.30 

-4.92 

4.92 

897.39 

-5.23 

5.23 

14 

71 .87 

1.137 

.692 

5 021 

98.18 

36.61 

36.61 

98.34 

36.84 

36.84 

15 

184.1 

1.264 

.6905 

30 400 

231.56 

25.78 

25.78 

228.88 

24.32 

24.32 

16 

684.5 

1.242 

.6927 

126 389 

716.47 

4.67 

4.67 

708.37 

3.49 

3.49 

17 

179.7 

1.307 

.6902 

31 608 

204.38 

13.73 

13.73 

202.27 

12.56 

12.56 

18 

659.6 

1.264 

.7188 

135 668 

591.70 

-10.29 

10.29 

688.99 

4.46 

4.46 

19 

71.86 

1.186 

.7158 

5 769 

77.12 

7.32 

7.32 

89.90 

25.10 

25.10 

20 

188.9 

1.324 

.7145 

33 921 

174.26 

-7.75 

7.75 

200.90 

6.35 

6.35 

21 

744.5 

1.244 

.6925 

166 931 

876.03 

17.67 

17.67 

863.18 

15.94 

15.94 

22 

240.2 

1.296 

.6901 

41 982 

262.82 

9.42 

9.42 

259.10 

7.87 

7.87 

23 

716.5 

1.238 

.6906 

124 455 

728.65 

1.70 

1.70 

710.73 

-.81 

.81 

24 

59.23 

1.203 

.6887 

5 190 

79.01 

33.40 

33.40 

77.94 

31 .59 

31.59 

25 

157.4 

1.335 

.688 

30 864 

184.48 

17.20 

17.20 

180.47 

14.66 

14.66 

26 

862.4 

1.193 

.6943 

115 964 

802.36 

-6.96 

6.96 

798.55 

-7.40 

7.40 

27 

157.1 

1.31 

.6897 

28 773 

189.14 

20.39 

20.39 

186.78 

18.89 

18.89 

28 

629.4 

1.242 

.6923 

161 542 

862.39 

37.02 

37.02 

848.76 

34.85 

34.85 

29 

237.15 

1.312 

.6898 

41 113 

244.85 

3.09 

3.09 

241.27 

1.59 

1.59 

30 

745 

1.249 

.6921 

143 053 

768.34 

3.13 

3.13 

756.45 

1.54 

1 .54 

31 

71.89 

1.205 

.6898 

5 999 

86.73 

20.65 

20.65 

86.07 

19.72 

19.72 

32 

200 

1.324 

.6896 

35 676 

211.41 

5.71 

5.71 

208.51 

4.25 

4.25 


ave. 


15.25 



15.26 


3 
























TABLE II 


TABLE III 


[Number of datasets = 32, Y = -2.78006 
-4.65772 XI -4.38901 X2 + 0.745065 
X3, Coefficiency of determination = 
0.974892] 


Y 

XI 

X2 

X3 

5.540596 

0.228489 

-0.37527 

10.09344 

4.913095 

.118937 

-.37062 

8.447628 

5.613638 

.249434 

-.37367 

10.23544 

5.575039 

.245922 

-.36976 

10.35866 

4.607467 

.150658 

-.37236 

8.735685 

5.859645 

.241690 

-.37367 

10.51680 

4.247780 

,161268 

-.37294 

8.287528 

5.293103 

.266739 

-.37498 

10.07382 

5.115175 

.106699 

-.37091 

8.709300 

5.696623 

.258124 

-.37367 

10.49174 

4.477223 

,155635 

-.37265 

8.569595 

5.542360 

.259976 

-.37411 

10.35427 

6.853193 

.158711 

-.36398 

11.70968 

4.274858 

.128393 

-.36816 

8.521384 

5.215479 

.234281 

-.37033 

10.32219 

6.528688 

,216722 

-.36715 

11.74711 

5.191288 

.267734 

-.37077 

10.36116 

6.491633 

.234281 

-.33017 

11.81796 

4.274719 

.170586 

-.33435 

8.660254 

5.241217 

,280657 

-.33617 

10.43178 

6.612712 

.218331 

-.36744 

12.02533 

5.481471 

.259282 

-.37091 

10.64499 

6.574378 

.213497 

-.37019 

11.73169 

4.081428 

.184818 

-.37294 

8.554488 

5.058790 

.288931 

-.37396 

10.33734 

6.759719 

.176471 

-.36485 

11.66103 

5.056882 

.270027 

-.37149 

10.26719 

6.444766 

.216722 

-.36773 

11.99252 


[Number of datasets = 32, Y = -1.10530 
-4.55257 XI + 0.736769 X2, 

Coefficiency of determination = 
0.973064] 


5.540596 

0.228489 

10.09344 

4.913095 

.118937 

8.447628 

5.613638 

.249434 

10.23544 

5.575039 

.245922 

10.35866 

4.607467 

.150658 

8.735685 

5.859645 

.241690 

10.51680 

4.247780 

.161268 

8.287528 

5.293103 

.266739 

10.07382 | 

5.115175 

.106699 

8.709300 

5.696623 

.258124 

10.49174 

4.477223 

.155635 

8.569595 

5.542360 

.259976 

10.35427 

6.853193 

.158711 

11.70968 

4.274858 

.128393 

8.521384 

5.215479 

.234281 

10.32219 

6.528688 

.216722 

j 11.74711 

5.191288 

.267734 

I 10.36116 

6.491633 

.234281 

11.81796 

4.274719 

.170586 

8.660254 

5.241217 

.280657 

10.43178 

6.612712 

.218331 

12.02533 

5.481471 1 

.259282 

10.64499 

6.574378 

.213497 

11.73169 

4.081428 

.184818 

8.554488 

5.058790 

.288931 

10.33734 

6.759719 

.176471 

11.66103 

5.056882 

.270027 

10.26719 

6.444766 

.216722 

11.99252 
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TABLE IV 


TABLE V 


[Nu = 13.865 (Tb/Tf )~-3.280 (Re)~0.3906 ...(18)] 



(18) 

Number 

Nu 

Tb/Tf 

Re 

Nucal . 

dNu(%) 

abs( ) 

1 


1.4364 

88 963 

362.74 

3.61 

3.61 

2 

391 .09 

1.4527 

102 357 

369.24 

-5.59 

5.59 

3 

397.01 

1.2892 

15 365 

260.41 

-34.41 

34.41 

4 

367.61 

1.4526 

115 255 

386.85 

5.24 

5.24 

5 

271.22 

1.2385 

7 853 

228.52 

-15.74 

15.74 

6 

435.62 

1.4362 

135 187 

427.35 


1 .90 

7 

384.29 

1.3643 

21 466 

246.45 

-35.87 

35.87 

8 

232.18 

1.3275 

12 616 


-5.66 

5.66 

9 

424.8 

1.4407 


357.59 

-15.82 

15.82 

10 

374.45 

1.4583 

129 718 

399.97 

6.82 

6.82 

11 

150.51 

1.3665 

12 563 

198.86 

32.12 

32.12 

12 

331.15 

1.4696 


369.61 

11.62 

11.62 

13 

103.59 

1.4221 

10 680 

163.75 



14 

305.2 

1.4615 

87 732 


11.68 

11.68 

15 

242.13 

1.4022 

13 899 

190.08 

-21.50 

21.50 

16 

196.34 

1 .2273 

8 630 

244.27 

24.41 

24.41 

17 

296.96 

1.5014 

91 088 

316.63 

6.62 

6.62 

18 

247.95 

1.4366 




25.38 

19 

177.72 

1.3432 

9 889 

191.64 

7.82 

7.82 

20 

338.56 

1.4676 



10.11 

10.1 1 

21 

260.23 

1.4272 

17 925 

198.12 

-23.87 

23.87 

22 

136.08 

1.3846 

10 878 


32.30 

32.30 

23 

404.85 

1.4496 

134 242 

413.40 

2.11 

2.1 1 

24 

258.12 

1.4363 


205.68 



25 

131.54 

1.4291 

4 915 

118.99 

-9.54 

9.54 

26 

330.2 

1.4543 

87 277 

345.70 

4.69 

4.69 

27 

260.23 

1.3892 

13 906 


-24.67 

24.67 

28 

171 

1.2539 

8 372 


31.58 

31.58 

29 

343.44 

1 .474 

128 536 

384.78 

12.04 



258.02 

1.4332 

19 375 


-21 .93 

21 .93 

31 

80.52 

1.48 

11 848 

149.60 

85.79 

85.79 

32 

311.63 

1.4806 

112 315 

359.72 

15.43 

15.43 

33 

287.5 

1.5903 

17 891 

138.82 

-51 .71 

51 .71 

34 

69.57 

1.4942 


136.90 

96.78 

96.78 


ave. 



22.73 


[Number of datasets = 34, 

Y = 2.629 -3.280 XI 
+0.3906 X2, coefficiency 

of determination = 0.788] 


Y 

XI 

X2 

5.858247 

0.362139 

11.39597 

5.968937 

.373423 

11.53622 

5.983961 

.254021 

9.639847 

5.907022 

.373355 

11.65490 

5.602930 

.213900 

8.968650 

6.076770 

.362000 

11.81441 

5.951397 

.310641 

9.974225 

5.447512 

.283297 

9.442721 

6.051618 

.365129 

11.38445 

5.925458 

.377271 

11.77311 

5.014029 

.312252 

9.438511 

5.802571 

.384990 

11.63590 

4.640440 

.352134 

9.276128 

5.720967 

.379463 

11.38204 

5.489474 

.338042 

9.539572 

5.279847 

.204816 

9.062999 

5.693597 

.406398 

11.41958 

5.513227 

.362279 

9.674074 

5.180209 

.295054 

9.199178 

5.824701 

.383628 

1 1 . 64646 

5.561565 

.355714 

9.793951 

4.913242 

.325411 

9.294497 

6.003516 

.371287 

11.80739 

5.553424 

.362070 

9.943140 

4.879310 

.357044 

8.500047 

5.799698 

.374524 

11.37684 

5.561565 

.328728 

9.540075 

5.141663 

.226258 

9.032648 

5.839012 

.387979 

11.76396 

5.553037 

.359909 

9.871738 

4.388505 

.392042 

9.379914 

5.741816 

.392447 

11.62906 

5.661222 

.463922 

9.792053 

4.242333 

.401590 

9.233079 
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TABLE VI 


19.402 (Tbr-0.538 (pTO.3593 (v)~0.9794 


Number 


116.1 

20.52 

39.22 

92.59 
41 

91.91 

21.69 

44.3 
102.87 

55.03 

96.7 

18.4 

37.23 

108.5 

35.27 

87.2 

21.56 

55.8 

102.8 

22.51 

47.11 

59 

66.63 

64.29 

32.65 

86.01 

22.49 

48.71 

51.99 

74.92 

28.69 
63.38 


1035.1 

1166.2 

1144.8 
1070 

1293.6 

1062.9 

1160.9 

1027.1 

917.8 

1243.9 
911 

1144.8 

1065.3 

917.4 
1213 

1066.2 

929.4 

1228.4 

1140.8 

1237.6 

1237.6 

940.2 

1242.6 

934.9 
1271 

855 

1303.3 

945.7 

1383.7 


85.99 

19.07 

65.23 

96.93 
69.33 

101.19 

20.07 

71.46 
105.41 

73.09 
94.89 

20.76 

72.96 

86.46 

68.97 

104.6 

22.06 

75.63 
101 .21 
21.6 

74.03 
61 .07 
67.65 
70.22 

21.3 

72.99 

19.2 
67.68 

18.95 

74.9 
20.75 

73.09 


37.56 

15.21 

23.53 

42.5 
26.62 

47.4 

19.2 

29.69 

53.96 

33.8 

42.47 

17.2 

26.6 

37.56 

23.53 

52.93 
21.84 
33.78 

48.03 
19.45 

30.08 

23.76 

27.43 

30.06 

22.22 
34.37 
15.31 
23.68 
22.26 

34.43 
19.35 
29.92 


7.7 

2.2 

5.6 


94.02 

49.43 
100.64 

24.08 

52.59 

120.69 

63.81 

93.63 

23.77 

48.59 

85.64 
44.16 

117.96 

30.71 

63.06 

105.94 

27.02 

55.47 

42.76 
48.87 
54.18 

30.43 
62.50 
20.42 

41.72 

30.77 

61.59 
26.25 
51.52 


dh(%) 

abs( ) 

-24. 

.40 

24, 

.40 

8. 

,33 

8, 

,33 

16. 

,47 

16. 

AT 

1 . 

,54 

1 . 

.54! 

20, 

,57 

20. 

.57 

9. 

.50 

9, 

.50 

11, 

.03 

11, 

.03 

18. 

,72 

18. 

,72, 

17, 

,32 

17. 

.32 

15. 

,96 

15. 

,96 

-3. 

.18 

3, 

.18 

29, 

.19 

29 

.19 

30, 

.52 

30 

.52 

-21 . 

.07 

21, 

.07 

25. 

.21 

25, 

.21 

35. 

,27 

35, 

.27 

42. 

.42 

42, 

.42 

13. 

.00 

13 

.00 

3. 

.06 

3 

.06 

20. 

,01 

20, 

.01 

17. 

,75 

17, 

.75 

-27. 

.53 

27, 

.53 

-26. 

.65 

26, 

.65 

-15. 

.73 

15, 

.73 

-6. 

.79 

6, 

.79 

-27. 

,33 

27, 

.33 

-9. 

,21 

9, 

.21 

-14. 

,35 

14, 

.35 

-40, 

,82 

40, 

.82 

-17. 

.79 

17, 

.79 

-8. 

.51 

8, 

.51 

-18. 

,71 

18. 

.71 













TABLE VII 


[h = 6.2315 ( Tb) ~-0 .450 (PTo. 0689 (v)~1.4327 ...(22)] 



(22) 

Number 

h 

Tb 

P 

V 

heal . 

dh(%) 

abs( ) 

1 

58.57 

1973.6 

170.93 

39.74 

57.03 

-2 

62 

2.62 

2 

50.59 

1681.6 

31.53 

31.72 

39.50 

-21 

93 

21.93 

3 

38.55 

1019 

16.72 

24.98 

33.64 

-12 

73 

12.73 

4 

20.18 

882.6 

17.25 

15.68 

18.46 

-8 

54 

8.54 

5 

63.61 

2294.1 

188.8 

43.79 

61.67 

-3 

05 

3.05 

6 

56.53 

1912.2 

37.17 

36.62 

46.32 

-18 

07 

18.07 

7 

40.63 

1341.1 

22.81 

28.84 

37.31 

-8 

16 

8.16 

8 

17.6 

971.7 

19.69 

18.1 

21.91 

24 

48 

24.48 

9 

68.59 

2095.4 

193.76 

47.99 

73.38 

6 

98 

6.98 

10 

58.62 

1875 

37.54 

40.14 

53.33 

-9 

02 

9.02 

11 

33 

1483.7 

25.99 

31 .61 

41.03 

24 

33 

24.33 

12 

20.25 

1056.8 

22.05 

19.84 

24.25 

19 

75 

19.75 

13 

81 .42 

2050.9 

196.21 

54.86 

89.82 

10 

32 

10.32 

14 

59.84 

1956.3 

40.52 

45.88 

63.70 

6 

45 

6.45 

15 

34.23 

1669.8 

30.26 

36.13 

47.61 

39 

09 

39.09 

16 

62.98 

1945.5 

168.34 

37.8 

53.38 

-15 

25 

15.25 

17 

53.43 

1629.1 

30.52 

31.61 

39.78 

-25 

56 

25.56 

18 

34.92 

1080.4 

17.71 

24.89 

32.73 

-6 

28 

6.28 

19 

71.34 

2179.8 

202.86 

54.96 

87.82 

23 

10 

23.10 

20 

59.15 

1929.1 

38.87 

45.96 

64.08 

8 

33 

8.33 

21 

64.27 

2168.2 

199.61 

47.77 

71.93 

11 

92 

11.92 

22 

26.15 

1601 

27.92 

31.46 

39.57 

51. 

,33 

51.33 

23 

64.81 

1853.5 

160.53 

37.94 

54.66 

-15. 

,66 

15.66 

24 

77.3 

2022.6 

181 .6 

43.79 

65.09 

-15. 

79 

15.79 

25 

82.98 

1541.8 

29.97 

36.62 

50.28 

-39. 

,41 

39.41 

26 

72.84 

2017.5 

186.56 

47.99 

74.44 

2. 

20 

2.20 

27 

86.25 

1515.5 

30.34 

40.14 

57.84 

-32. 

94 

32.94 

28 

85.94 

1984 

189.81 

54.86 

90.96 

5. 

85 

5.85 

29 

80.86 

1647.3 

34.12 

45.88 

68.01 

-15. 

89 

15.89 

30 

52.84 

1316.6 

23.86 

36.13 

52.12 

-1 . 

35 

1.35 

31 

66.53 

1880.8 

162.74 

37.8 

54.07 

-18. 

73 

18.73 

32 

76 

2093.8 

194.86 

54.96 

89.18 

17. 

34 

17.34 

33 

35.65 

1420.9 

25.05 

36.2 

50.68 

42. 

15 

42.15 

34 

67.17 

2107.4 

194.01 

47.77 

72.71 

8. 

25 

8.25 

35 

26.15 

1601 

27.92 

31.46 

39.57 

51. 

33 

51.33 


ave 



17.83 
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FIGURE 1. - INSTRUMENTED OMC ENGINE 
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FIGURE 2. - COOLANT FLOW THROUGH ENGINE 
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FIGURE 3. - TIME AVERAGE WALL TEMPERATURE. GAS SIDE WALL TEMPERATURE VERSUS 
C.A. 
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FIGURE - STRAIN GAGES AND THERMOCOUPLES LOCATION. 
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FIGURE 5. - AVERAGE HEAT FLUX PER UNIT AREA. HEAT FLUX PER UNIT AREA VERSUS 
C.A. 



FIGURE 6. - PRESSURE TRANSDUCER LOCATIONS. 
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FIGURE 7. - A-D REAL TIME DATA RECORDING SYSTEM 




NASA 

National Aeronautics and 
Space Administration 

1. Report No. 

NASA TM- 1001 24 

4. Title and Subtitle 


Report Documentation Page 


2. Government Accession No. 


3. Recipient’s Catalog No. 


5. Report Date 


Regressed Relations for Forced Convection Heat 
Transfer in a Direct Injection Stratified Charge 
Rotary Engine 


7. Author(s) 


Chi M. Lee and Harold J. Schock 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


6. Performing Organization Code 


8. Performing Organization Report No. 


E— 3689 

10. Work Unit No. 

505-62-11 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 

14. Sponsoring Agency Code 


15. Supplementary Notes 

Prepared for the 1988 International Congress and Exposition sponsored by the 
Society of Automotive Engineers, Inc., Detroit, Michigan, February 29 - March 4, 
1988. Chi M. Lee, NASA Lewis Research Center; Harold J. Schock, Michigan State 
University, East Lansing, Michigan 48824. 

16. Abstract 

Currently, the heat transfer equation used in the rotary combustion engine (RCE) 
simulation model is taken from piston engine studies. These relations have been 
empirically developed by the experimental input coming from piston engines whose 
geometry differs considerably from that of the RCE. The objective of this work 
was to derive equations to estimate heat transfer coefficients in the combustion 
chamber of a RCE. This was accomplished by making detailed temperature and pres- 
sure measurements in a direct injection stratified charge (DISC) RCE under a 
range of conditions. For each specific measurement point, the local gas velocity 
was assumed equal to the local rotor tip speed. Local physical properties of the 
fluids were then calculated. Two types of correlation equations were derived and 
are described in this paper. The first correlation expresses the Nusselt number 
as a function of the Prandtl number, Reynolds number, and characteristic tempera- 
ture ratio; the second correlation expresses the forced convection heat transfer 
coefficient as a function of fluid temperature, pressure, and velocity. 


17. Key Words (Suggested by Author(s)) 

Direct disk rotary engine; 

Forced convection heat transfer 
coefficient; Multiple regression 
method 


18. Distribution Statement 

Unclassified - Unlimited 
Subject Category 07 


19. Security Classif. (of this report) 

Unclassified 


20. Security Classif. (of this page) 

Unclassified 


21. No of pages 


NASA FORM 1626 OCT 86 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 





